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^ | I Abstract. Heterodyne, eight-port homodyne and six-port homodyne detectors 

^vq ' belong to the class of two-photocurrent devices. Their full equivalence in probing 

CNI ■ radiation field has been proved both for ideal and not fully efficient photodetectors. 

The output probability distribution has been also evaluated for a generic probe mode. 
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1. Introduction 

vo . , 

Q\ ' In order to gain information about a quantum state of the radiation field one has to 
measure some observable. The measurement process unavoidably involves some kind of 
interaction, thus coupling the mode under examination to one or more other modes of 



a: 



the field. One has to admit that, in general, the measured observable is not defined on 
the sole Hilbert space of the signal mode On the contrary, it reflects properties of 
the global state, possibly entangled, coming from the interaction among the signal mode 
and the set of the probe modes. Sometimes, it is possible to get rid of these probe modes, 
^ ■ so that the statistics of the output is related only to the quantum statistics of the signal 
mode. This, as an example, is the case of balanced homodyne detectors @]. In that 
case, in fact, an appropriate rescaling of the output photocurrent allows to completely 
neglect the local oscillator in the definition of the measured observable, which simply 
results to be proportional to the signal field quadrature a(0) = ^(a^e 1 ^ + ae" 1 ^). 

However, in more general cases this procedure cannot be pursued and the output 
statistics always remind us the way we are probing the signal under examination. This is 
a common feature of three relevant detection schemes in quantum optics, which provide 
different setup for jointly measuring a couple of photocurrents. They are the heterodyne 
0, f|, ||, the eight-port homodyne H 0, ||, || and the recently introduced six-port 
homodyne detectors ]ID|, |TT| . 

It is a purpose of this paper to show that, although they involve very different ways 
of coupling signal to probe modes, they provide the same information on the signal under 
examination. Actually, we will prove their full equivalence by demonstrating that their 
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output photocurrents have the same operatorial form. More surprisingly, this remains 
true also when taking into account the inefficiency of photodetection, even though the 
latter occurs at very different stages in the three detection schemes. 

The paper will be organized as follows. In the next section we briefly review how 
to describe inefficient photodetection in terms of beam splitters and ideal detectors. In 
the three sections 0|J and |5] we will examine the three detection schemes respectively, 
in order to show that all of them jointly measure the real and the imaginary part of the 
complex photocurrent [T2j 

-§ = + 6+, (1) 

being a the signal mode and b a probe mode. The inefficiency of realistic photodetectors 
will be also taken into account, in order to show it does not affect the equivalence of 
the considered detection schemes. 

In section refs:pom we will derive the output statistics of considered schemes as 
members of the more general class of two-photocurrent devices. Section |7] will close the 
paper with some concluding remarks. 



2. Inefficient photodetection 

The final stage of any detection scheme is represented by photodetection, namely 
counting photon through their conversion to electric pulses. Let us consider a light beam 
p entering a phototube which converts to electric pulses a fraction ( of the incoming 
photons. By keeping open the detector window for a time interval T, the probability 
P m (T) of counting m photons is given by (T^j 



P m (T) = Trip p™™ exp[-CJ(T)r]:| , (2) 



/(T) = ^/ T E H (r,t).E (+) (r,^. (3) 



where : : denotes normal ordering of operator and I(T) is the beam intensity 

1 Jo 

E*' ^ (r, t) denotes the positive (negative) frequency part of the field. For a single-mode 
field excited in a stationary state Eq. |] can be rewritten as 

C = Tr|p^^exp(-Wa)| , (4) 

being [a, a^] = 1 the single mode field operator and r] = (chu/V the global quantum 
efficiency of the photodetectors. For unit quantum efficiency this coincides with the 
actual photon number distribution of the state under examination 

Pm = Pmm = (m\p\m) , (5) 
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whereas, in the realistic case of non-unit quantum efficiency Eq. (Hf) becomes a binomial 



convolution In formula 

Pl = flPnni" 1 ) V m (l-V) n ~ m - (6) 



HI 

n=m \ 



Let us now consider the scheme in Fig. |T|. The signal mode a is impinged in a beam 
splitter of transmissivity r whose second port b is placed in the vacuum. A perfect 
(77 = 1) photodetection on the exiting mode reveals m photons with a probability 

P m = Tr(y r (p®\0)(Q\)U}\m)(m\®i) , (7) 

being 1 the identity operator on the second input of the beam splitter and 



tl T = exp I i arctan 4 / (cftb — atf) 



the evolution operator of the beam splitter. A straightforward calculation shows that 

P m =jrPnn( U )(l- r) n - m T m . (8) 



n—rn 



Eq. (H) coincides with Eq. (^) for r = 77. This means that a photodetection process with 
efficiency rj is equivalent to a perfect photodetection process performed on the original 
signal mixed with vacuum by a beam splitter with a value of the transmissivity equal to 
the quantum efficiency |15|, |l6f. In the following we will adopt this equivalent scheme. 



3. Eight-port homodyne detector 

Eight-port homodyne detector is known for a long time for the joint determination of 
phase and amplitude of the field in the microwave domain. It was introduced in the 
optical domain by Walker and Carrol || and successively analyzed by different authors 
0, 1, 1, 0, 0, |Tj, 0, H, H, H H, H M ■ A schematic diagram of the experimental 



setup is reported in Fig. [|. There are four balanced beam splitters whereas a 7r/2 phase 
shifter is inserted in one arm. The four input modes are denoted by a^, k = 1,...,4 
whereas the detected output modes are denoted by bk, k = 1,...,4. There are four 
identical photodetectors whose quantum efficiency is given by rj. The noise modes used 
to take into account inefficiency, according to the scheme of the previous section, are 
denoted by u iy % = 1, ...,4. We consider a>i as the signal mode, whereas a 2 is referred to 
be the idler of the device. The mode 03 is unexcited, whereas a 4 is placed in a highly 
excited coherent state \z) provided by an intense laser beam (local oscillator). The 
detected photocurrents are -4 = b\,bk, k = 1, ...,4 which form the eight-port homodyne 
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observables 



Z, 



h-h 
2t]\z\ 

h-h 
2t]\z\ 



(9) 



The latter are derived by rescaling the difference photocurrent, each of them obtained 
in an homodyne scheme. For this reason eight-port homodyne is known also as double- 
homodyne detection. In Eq. (ffl) rj denotes the quantum efficiency of the photodetectors 
whereas \z\ is the intensity of the local oscillator. In order to obtain Z\ and Z 2 in terms of 
the input modes we first note that the input-output mode transformation is necessarily 
linear, as only passive components are involved in the detection scheme |28|, |29j . Thus, 
we can write 

4 

b k = ]T M kl a t , (10) 
1=1 

where the transformation matrix can be computed starting from the corresponding 



transformations for the beam splitters and the phase shifter |p2 

1111 

1 i -1 -i 
1 -1 i -1 
I -i -I i 



M 



1 

74 



Eqs. ([H]) and (|TT|) together with the equivalent scheme for the inefficient detector leads 
to the following expression for the output modes 
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h = yfrj [at + a 2 + a 3 + a 4 ] + yj 1 
sfr\ [ax + ia 2 — a 3 — ia^ 



V u 2 



a 2 + ia 3 — a 4 J + \J1 — 77 w 3 
ia 2 — 03 + ia^} - 



l-rju 4 . (12) 
and considering the limit of highly excited local 



Upon inserting Eqs. ( [12] ) in Eq. 1 
oscillator we obtain the eight-port photocurrents in terms of the input modes 



Z x = ai(0) + a 2 (0) 



1 — 77 



«i(0)-«2(0)] + O[ 



Zo 



-ai(vr/2) + a 2 (n/2) + 



^[n 3 (7r/2)-n 4 (7r/2)] + 0[^]. 

Tj \Z\ 



(13) 



In Eq. ( p!3|) a(0) denotes a quadrature of the field. The complex photocurrent 
Z = Z\ + iZ 2 is given by 



Z = ax + a 2 , 



(14) 
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for unit quantum efficiency, whereas for non unit quantum efficiency it becomes a 
Gaussian convolution of Eq. (fH]), we will consider this point in detail in ^. 

It is worth noticing here that the mode transformation defined by Eqs. ( [T0| ) and ([TT|) 
is distinctive for a canonical 4 x 4-port linear coupler as defined in Refs. |31| . It has 
been rigorously shown that a N x iV-port linear coupler can always be realized in 
terms of a number of beam splitters and phase-shifters. However, this implementation 
is, in general, not unique. The interest of eight-port homodyne scheme lies in the fact 
it provides the minimal scheme for realizing a 4 x 4-multiport |B3. If one considers the 



multiport as a given black-box device the eight-port homodyne scheme can be depicted 
as in Fig. |3|. This will facilitate the comparison with the six-port homodyne detection 
presented in Section ||. 

4. Heterodyne detector 

Heterodyne detection scheme is known for a long time in radiophysics. It has been 
introduced in the domain of optics || £|, ||, |34| in order to describe the joint measurement 



of two conjugated quadratures of the field. The term 'heterodyne' is used as the involved 
field modes are excited on different frequencies. 

In Fig. H we show a schematic diagram of the detector. We denote by E5 the signal 
field, whereas E/x> describes the local oscillator. The field E £ accounts for the losses 
due to inefficient photodetection. The input signal is excited in a single mode (say a) 
at the frequency u, whereas also the local oscillator is excited in only one mode at the 
frequency uj . This local oscillator mode is placed in a strong coherent state \z) by means 
of an intense laser beam. The beam splitter has a transmissivity given by r, whereas 
the photodetectors shows quantum efficiency rj. The heterodyne output photocurrents 
are given by the real Z\ and the imaginary Z2 part of the complex photocurrent Z. The 
latter is obtained after the rescaling of the output photocurrent /, which is measured 
at the intermediate frequency u>i = uj — uj . By Fourier transform of Eq. (||]) we have 

= [ dJ E^~V + uj) Eq V) , (15) 

JR. 

being Eq ^ positive and negative part of the output field. In terms of the input fields 
Eq. ( [T3|) can be written as 

I{uj) = fdu/ U^s~V + wj) + V^(1-t)e1oV + wj) + 0^?Ei~V + u T ) 

J Ft L 



(16) 



Heterodyne photocurrent is obtained by the following rescaling 



Z = lim = with — r constant . (17) 

r-»i,[z|— 00 |z|^J T (l - r) 
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Physically this definition corresponds to require a very intense local oscillator, which 
however is allowed only for a little mixing with the signal mode |53]. In this limit only 
terms containing the local oscillator field E^q(o;o) at the frequency ujq can survive in 
Eq. fll6|), so that we have 



Z = Z\ + iZ 



2 • 



;is) 



where 



Zx = a(0) + c(0) + 



(ui(0) - u 2 (0)) 



2o 



a(vr/2) + c(?r/2) + 



In writing Eq. fll9"|) we have substituted 



c <— E^(2o;o — c<j) 



'1 — rj 



(u 1 (n/2)-u 2 (ir/2)) . 



(19) 



It! 



,(+) 



KM 



(20) 



u 2 <- E^ +) (2^ - oS) 

for the relevant modes involved. Provided u\ and u 2 to be noise modes placed in the 
vacuum the expression ( jl9|) for the heterodyne photocurrents leads to the identical 
output statistics of eight-port homodyne photocurrents, the mode c playing the role of 
the idler of the device. The full equivalence of the two detection schemes has been thus 
proved. 



5. Six-port homodyne detector 



A linear, symmetric three-port optical coupler is a straightforward generalization of 
the customary lossless symmetric beam splitter. The three input modes a«, i = 1,2,3 
are combined to form 3 output modes bj, j = 1,2,3. In analogy to lossless beam 
splitters, which are described by unitary 2x2 matrices |36|, any lossless triple coupler 



is characterized by a unitary 3x3 matrix |37|, |38| . For the symmetric case we have the 
form 

/ll 1 \ 



T = -^= 



1 exp{if} exp{-zf } 



V 



(21) 



/ 



1 exp{-if } exp{if } 
where each matrix element 7^ represents the transmission amplitude from the z-th input 



port to the j-th output port, that is bj 



1 TjkCLk- 
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Such devices have already been implemented in single-mode optical fiber technology 
and commercial triple coupler have been available for some time [p9| . Any triple 
coupler can be also implemented by discrete optical components using symmetric beam 
splitters and phase shifters only |37|. As it has already mentioned in Section this 
is due to remarkable mathematical fact that that any unitary M-dimensional matrix 
can be factorized into a sequence of 2-dimensional transformations plus phase-shift 
Moreover, this decomposition is not, in general, unique. In Fig. [5] a possible 



implementation of a triple coupler is schematically reported. Experimental realizations 
of triple couplers has been reported for both cases, the passive elements case and the 



optical fiber one |7|, H, g(J ||. 

Let us now consider the measurement scheme of Fig. || The three input modes are 
mixed by a triple coupler and the resulting output modes are subsequently detected 
by three identical photodetectors. The measured photocurrents are proportional to J n , 
n — 1, 2, 3 given by 

L = b% = \ £ exp {%9 n {l - k)} a\ ai , e n = ^-(n-l). (22) 
6 k,i=i 6 

After photodetection a Fourier transform (FT) on the photocurrents is performed 

1 3 . 

l s =FT(/ 1; I 2 ,/3) = _£/ n exp{-i0 n (s-l)} , S = l,2,3. (23) 

V 3 n= i 

This procedure is a straightforward generalization of the customary two-mode balanced 
homodyning technique. In that case, in fact, the sum and the difference of the two 
output photocurrents are considered, which actually represent the Fourier transform in 
a two-dimensional space. By means of the identity 



S 3 (s - 1) 



lX>p{4pi( S -l)} , (24) 



for the periodic (modulus 3) Kronecker delta 63 , we obtain our final expressions for the 
Fourier transformed photocurrents 

It = -j= \a\ai + a\a 2 + a\a^ , (25) 
\a\a 2 + ajos + 4ai} , (26) 



I, 
X, 



^= \a\az + a|ai + 4a 2 } . (27) 



X\ gives no relevant information as it is insensitive to the phase of the signal field, whereas 
X2 and X3 are hermitian conjugates of each other and contain the relevant information 
in their real and imaginary part. In the following let us assume a\ is the signal mode 
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and a 2 is fed by a highly excited coherent state \z) representing the local oscillator. For 
large \z\ the output photocurrents are intense enough to be easily detected. They can 
be combined to give the reduced photocurrents 

ii = v^^p = ai(o) + a 3 (o) + o[h 

Z 2 = v^^p = -a 1 (7T/2) + a 3 (vr/2) + 0[h , (28) 

which we refer to as the triple homodyne photocurrents. Again the complex photocurrent 
Z = Z\ + iZ 2 has the form Z = a\ + a\, being a\ the signal mode and 03 the idler of 
the device. 

When accounting for the non unit quantum efficiency rj of the photodetectors the 
output modes are written as 



bj = y/rj T jk a^j 



-rjUj j = 1,2,3 
so that the reduced photocurrents are now given by 



Zi = = + a 3(o) + J— Mo) - 63(0)] + o[h 

2i]\z\ y 7] \z\ 

Z 2 = v^^-jjp = -ai(vr/2) + + J— ^ [ui(tt/2) - u 3 (tt/2)] + 0[— } . (29) 

When, as it is the case, the modes Uj are placed in the vacuum the six-port photocurrents 
in Eq. ( |2"9"D leads to the same statistics of the eight-port photocurrents in Eq. fllBp. 
Indeed, they describe different devices leading to the same amount of information on 
the signal mode a\. Some comments are, however, in order. By comparison of Fig. [3] 
and Fig. || it appears obvious that six-port homodyne is an optimized scheme relative to 
the eight-port one. One mode less is needed to reach the same amount of information, 
thus decreasing the possible sources of noise. The reason for this lies in the final stage 
of the two schemes. The Fourier transform of the six-port photocurrents, in fact, is a 
more effective procedure relative to the double-homodyne final stage of the eight-port 
scheme. This is related to the noise suppression mechanism of homodyne detection. The 
latter, in fact, is generalized to the multi-mode case by the Fourier transform rather than 
duplication of the original two-mode scheme. 



6. Output statistics from a two-photocurrent device 



In this section we analyze with some detail the output statistics of an abstract two- 
photocurrent device. The latter is characterized by the joint measurement of the real 
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Z\ and the imaginary Z 2 part of the complex photocurrent 

Z = a + 6 f , (30) 

when equipped with perfect photodetectors. On the other hand, in the realistic case of 
inefficient photodetection the photocurrent is given by 



Z = a + tf + J ^(ut + ul) . (31) 




In Eqs. (|30|) and (^l|) a and b denote single mode radiation field which can be excited 
in any quantum state. The operators U\ and u 2 also denote single mode radiation field, 
however strictly placed in the vacuum states. They are used to simulate losses due 
inefficient photodetection, according to the equivalent detection scheme of Section 0. In 
the following we will refer to the mode a as the signal mode which is under examination, 
whereas the mode b is in a known and fixed state, playing the role of the probe mode of 
the device. This is only for sake of convenience. It is obvious that the roles of the two 
modes can be interchanged and any argument can be reversed considering the mode b 
as a signal mode. According to this scheme we denote the outcome probability density 
distribution by K b (a,a). The latter describes, in the complex plane, the state of the 
mode a as probed by the mode b. Indeed, different choices for the probe mode lead to 
very different features in the probability distribution. 

Each experimental random outcome (z\, z 2 ) from the joint measurement of Z\ and 
Z 2 can be considered as a point z in the complex plane. On the other hand, the output 
photocurrent Z is expressed as a sum of different contributions coming from different 
modes. Therefore, it appears intuitively rather obvious that the resulting probability 
distribution will be given by a convolution. To be more specific let us start by considering 
the ideal case of unit quantum efficiency 77 = 1. We write the probability distribution 
Kb(a, a) as the Fourier transform 

K b (a,a)= I ^e^ 5 5(7,7), (32) 

of the characteristic function 

5(7, 7) = Tr {f3 exp [7 Z - >yZ^ } , (33) 

being q the global density matrix describing both modes a and b. We consider probe 
mode to be independent on the signal mode, so that the input mode is factorized as 

Q = Qa® Qb- (34) 

Upon substituting Eqs. ([3(]) and ( |3"1| ) in Eq. ( |3"3"D we are able to write the characteristic 
function 5(7, 7) product 

5(7,7) =T±{Qa®Qb A*(7)®A,(-7)} = Xa(7,7)Xfc(-7,-7), (35) 
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being 

D(j) = {pexp ^a) — 7a 
the displacement operator and 

x{l,l) = Tr {^(7)} 

the single mode characteristics function, the latter entering in the definition of the 
Wigner function of a single mode radiation field W\ f43l E4 



W(a,a) = J x(A, A) exp |Aq; — Xaj . (36) 



We now insert Eq. (pop in Eq. (|3^). By means of Eq. (p6|) and using the convolution 
theorem we arrive at the final result 

K b (a, a) = W a (a, a) * W b (-a, -a) = [ ^ W 6 (a + /3, a + 0) W a (A P) , (37) 



the symbol * denoting convolution. Eq. ( p7|) provides a justification for referring to 
the mode b as the probe of the device. In fact, it acts as a filter on the signal' Wigner 
function. The latter is not a genuine probability distribution, as it can be negative when 
describing quantum interference. Thus, it cannot be directly sampled by experiments. 
However, the convolution (|37| ) make it more regular, leading to K b (a,a) which is a 
measurable distribution. 

Phase space density as in Eq. fl3"7p have been already introduced by Wodkiewicz 
45|, f46l |47| to account for the effect of the measuring apparatus in a joint measurement 



of conjugated variables. Originally, they have been termed phase-space propensities. 
More recently, they also have been used in entropic description of quantum mechanical 
states [^, Two-photocurrent devices appear as the natural setup to start from, in 
order to experimentally access such kind of phase-space distribution. 

Wigner function, though it contains a complete description of the quantum state, 
cannot be directly measured. The phase distributions K b (a, a) are smoothed version 
of it, corresponding to the occurrence of additional noise of purely quantum origin. 
Indeed, two-photocurrent devices provide the generalized joint measurement of position 
and momentum, thus unavoidably introducing additional noise by first principles 
f34| , |50| , [51]. The crucial point to notice here is that they offer the remarkable possibility 
to manipulate this quantum noise. As it emerges from Eq. (|37|) it can be redirected 
in the desidered region of the complex plane by suitable choice of the probe mode, 
according to the kind of information of which is of interest. 

In order to incorporate the effects of inefficient detection we first note that Eq. (pl|) 
differs from Eq. ( pCf ) by two additional additive terms. Therefore, further convolutions 
to the original Wigner function are expected. Indeed, the characteristic function of the 
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whole device is now expressed as 

2^(7> l) = Xa(l, l) Xb(-1, -7) x 

x , , 

77 y 77 y \ 1 t) V 7 ? 
where characteristic function for the noise modes can be easily evaluated as 



(3? 





exp 



l-^i ,2 

7 



277 



J = 1,2 



(39) 



Upon inserting Eqs. (|3^) and (|39|) in Eq. fl32|) we obtain our final result, namely the 
output probability distribution K bri (a,a) of a two-photocurrent device equipped with 
photodetectors of quantum efficiency 77 



c K 



K b ((3,(3) exp 



1 — 77 



|a-/3|' 



(40) 



being ^1(0;, a) the probability distribution obtained for ideal photodetection and 

G v (a, a) = exp 



V 1 ,2 
a 



1 — 77 

the filter function summarizing the effects of inefficient photodetection. 
7. Conclusions 

In this paper we have introduced the class of two-photocurrent devices. This kind 
of detectors are characterized by the fact that they jointly measure the real and the 
imaginary part of the complex photocurrent Z = a+tf , where a and b describe two single 
mode of the radiation field. Eight-port homodyne, heterodyne and six-port homodyne 
detectors belong to this general class, thus their fully equivalence in probing radiation 
field has been proved. It has also been proved that this equivalence still holds when the 
inefficiency of the photodetection process is taken into account. This is an interesting 
and unexpected result, as photodetection takes place at very different stages in the three 
detection schemes. 

The three schemes analyzed in this paper are equivalent from the point of view of 
provided information on the measured signal. Nevertheless, they have different physical 
implementations which have to be compared. We pointed out the advantages of the 
six-port homodyne scheme in comparison with the customary eight-port one. Actually, 
it provides the minimal scheme to access generalized phase space distribution [E2| . 

In a two-photocurrent device a generalized joint measurement of position and 
momentum is performed on the signal mode. This results in a smoothing of the signal 
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Wigner function to a measurable distribution, which represents the output probability 
distribution of the measurement. Some additional noise is unavoidably introduced, 
according to the Heisenberg principle for joint measurement. However, the filtering 
process has been shown to be a convolution with the Wigner function of the probe 
mode. Therefore, it is possible to manipulate and redirect the noise. A suitable choice 
of the probe mode enhances different features of the signal' phase space distribution, 
according to the kind of desidered information. 
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Figure 2. Schematic diagram of an eight-port homodyne detector. 
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Figure 4. Schematic diagram of a heterodyne detection. Relevant modes are explicitly 
pointed out. 



17 




Figure 5. Realization of a triple coupler in terms of 50:50 beam splitters (BS) and 
phase shifters '<p\ In order to obtain a symmetric coupler the following values has to 
be chosen: ipi — arccos(l/3) and = <pi/2. 




Figure 6. Outline of triple coupler homodyne detectors: The hexagonal box 
symbolizes the electronically performed Fourier transform (FT). 



